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We consider predictions for structure formation from modifications to general relativity in which 
the Einstein-Hilbert action is replaced by a general function of the Ricci scalar. We work without 
fixing a gauge, as well as in explicit popular coordinate choices, appropriate for the modification of 
existing cosmological code. We present the framework in a comprehensive and practical form that 
can be directly compared to standard perturbation analyses. 

By considering the full evolution equations, we resolve perceived instabilities previously suggested, 
, and instead find a suppression of perturbations. This result presents significant challenges for 

■ agreement with current cosmological structure formation observations. 

' The findings apply to a broad range of forms of f{R) for which the modification becomes important 

^Nj , at low curvatures, disfavoring them in comparison with the ACDM scenario. As such, these results 

t I . provide a powerful method to rule out a wide class of modified gravity models aimed at providing 

' an alternative explanation to the dark energy problem. 

0> ' I- INTRODUCTION 

' Approaches to resolving the late-time acceleration of the universe may be divided into three broad classes (see, 
^ , for example, [1] for a brief review). Perhaps the simplest possibility is that there is some, as yet undiscovered, 
^-H ' consequence of our existing model of gravity and matter that leads to acceleration during the current epoch. Included 
PsJ , in this category is the existence of a tiny cosmological constant and the possibility that the backreaction of cosmological 
' perturbations might cause self-acceleration. A second option is the idea that a new dynamical component exists in 
] cosmic energy budget. Included here are sources of energy density modeled by a scalar field, usually referred to as 
dark energy. Of equal interest, however, is that General Relativity (GR) can be modified in the low curvature regime, 
to admit self-accelerating solutions in the presence of negligible matter [2-21]. 

We consider, in this paper, the class of modified gravity models in which the gravitational action contains a general 
' function f{R) of the Ricci scalar. For such models, the analysis of the background cosmological evolution can be 
i^- significantly simplified by performing a conformal transformation on the metric. Such a transformation maps from 
Q ' a frame in which the gravitational action and resulting field equations are modified from GR, commonly called the 
^ ' Jordan frame, to a frame in which the gravitational action for the new metric is the Einstein-Hilbert one, commonly 
^2 ' called the Einstein frame. In this new frame, the matter fields couple non-minimally to the new metric (matter no 
longer falls along geodesies of the new metric) and an extra degree of freedom now manifests as a new scalar field. 
These couplings affect the evolution of perturbations in a potentially observable way [22]. 

This paper focuses on structure formation in general f{R) theories and comparison of the predictions with observa- 
tions. This complements previous work investigating the background evolution [23-26] and solar system implications 

■ [26-33] for f{R) modifications. As part of this analysis we derive the linear cosmological perturbation theory for 
a general form of f{R) considering only scalar perturbations, since the tensor and vector modes are unaffected by 
f{R) modifications to gravity. We present our equations without fixing a gauge, as well as in two explicit coordinate 
choices: the conformal Newtonian and synchronous gauges. These equations are presented in a comprehensive and 
practical form that can be directly compared to standard perturbation analyses (e.g. [34, 35]) and are applicable to 
the modification of existing cosmological codes such as C AMB [36] . This work complements alternative formulations 
in the Palatini approach [37], for the metric approach in the Einstein frame [38] and, more recently, in terms of 
"frame- independent" variables [39]. 

In section II we give an overview of the f{R) action and in section III discuss the conformal transformation used 
to express f{R) theories in a frame with Einstein field equations with non-minimally coupled matter. In sections IV 
and V we present the main analytical results of the paper. In section IV we present full perturbation equations for 
f{R) theories in the Jordan frame. In section V we examine the general behavior of late time structure formation in 
both Einstein and Jordan frames specifically for f{R) theories that allow late-time cosmological acceleration solutions. 
This leads to important constraints on a large class of f{R) theories dominating in the low curvature regime. For 
concreteness we present the predictions for two specific examples. We also outline the origin of the apparent zeroing of 
density perturbations when approximating the evolution equations in the Jordan frame [40] , and resolve the matter by 
considering the full Einstein frame equations. Finally in section VI we summarize our findings and discuss implications. 
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II. DESCRIPTION OF f{R) GRAVITY MODELS 

The class of theories on which we focus has action in the Jordan frame 

.5 = d'xv^ [R + f{R)] + J d\^C^[xi,g^u] , (1) 

where = SttG and the function /(i?) is a general function of the Ricci scalar, R. The matter Lagrangian, Cm, 
is minimally coupled and therefore the matter fields, Xi) fall along geodesies of the metric g^i,. The field equations 
obtained from varying the non- minimally coupled gravity action (1) with respect to g^i, are 

(1 + /fl) R^^u - Ig^^u {R + f ) + - V^V.) Jr = , (2) 

where we have defined fa = df/dR. We assume matter to behave as a perfect fluid, with energy- momentum tensor 

T^. = {p + P)U^U,+pg^^ , (3) 

where is the fluid rest-frame four-velocity, p is the energy density and P is the pressure. We relate P and p via 
p = wp, where w is the equation of state parameter (for pressureless matter w = and radiation to w = 1/3). 

When considering the background cosmological evolution, we take the metric to be of the flat Robertson- Walker 

form, ds'^ = a'^{T){—dT'^ + dx^), with a(T) the scale factor and r conformal time, in terms of which the curvature 
scalar satisfies R = 6a" /a^. Here a prime denotes differentiation with respect to conformal time, t, H = a' /a is the 
equivalent Hubble expansion rate. 

In the Jordan frame, the metric is that of standard GR minimally coupled to matter. Hence the stress tensor and 
its conservation laws will remain the ones of standard GR. The continuity equation is the usual 

p' + m{p + P) = {) . (4) 

However, since we have modified the gravitational action, there are extra terms appearing in the Einstein equations. 
In particular, for our cosmological ansatz, the Priedmann equation becomes 

(1 + + ^/ - ^fn + HfR = ^a^p (5) 



and the acceleration equation is 



^-{l + f^)H' + a'L+Hf'R + l-f'A = -^a\p + W). (6) 



III. MAPPING TO THE EINSTEIN FRAME 

Although the background cosmological behavior can be chosen by an appropriate choice of f{R), this freedom 
comes at the expense of unfamiliar dynamical equations. There exists, however, a complementary, and sometimes 
conceptually simpler, way in which to approach /(i?) modifications to GR. It is possible to perform a conformal 
transformation on the metric so as to render the gravitational action in the usual Einstein Hilbert form of GR. The 
price one pays for this simplification is a non- minimal coupling between matter fields and the new metric [41-43], 
as well as the appearance of a new scalar degree of freedom evolving under a potential determined precisely by the 
original form of the /(i?) coupling in the Jordan action. 

Using the approaches of Chiba [44] and of Magnano & Sokolowski [46], following [41], we recast the gravitational 
action (1) into a dynamically equivalent form by introducing an intermediate scalar field The equivalent action is 
[46] 

^ = ^ Jd'xV^[{^ + f{^)) + {l + U){R-n 

d'^Xy/^Cm[Xi,gfJ.>^] , (7) 



where /$ = 9//c5$. One can verify that, if d^f/d^^ ^ 0, the field equation for $ is = which reduces (7) to the 
original action. Next consider the conformal transformation 

~g^, = e2-(-")5^. , (8) 
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such that the function a; (a;") satisfies 



e-'"(l + /ii) = l . 



(9) 



With this choice of lj the action (7) transforms into an action with the usual Hilbert-Einstein form for gravity. If 
we now define the scalar field (j) = 2uj/pK, where (3 = -^2/3, the resulting action becomes 



1 

2^ 



d x\J —g R - 



-^r''{^,.4>)y.4>-v{4>) 



where the potential V{(j}) is determined entirely by the original form (1) of the action and is given by 

2k^ (1 + fny 

The Einstein-frame line element can be written in familiar Fricdmann-Robertson- Walker (FRW) form as 

dS"^ =d'^{-dT^ +dx'^) , 



(10) 



(11) 



(12) 



where the Jordan and Einstein metrics arc related through the conformal transformation d^ = e^'^'^a?. It is also 
convenient to define an Einstein-frame matter energy-momentum tensor by 



where % = e^'">'/'^ U^, p = e'^'^'^'I'p and P = e-^'^'"t'P. 

The equations of motion obtained by varying the action with respect to the metric g^v are 



(13) 



(14) 



and are more familiar than those in the Jordan frame, although there are some crucial distinctions. Most notably, 
in this frame test matter particles do not freely fall along geodesies of the metric g,^v, since the scalar field is also 
coupled to matter. 

The remaining equations of motion, for the scalar field and for the perfect fluid matter, are given respectively by 



f + 2n(t)' + a'V^ = ^KPd\p - 3P) , 

p' + snip + P) = -\i^P4>\p - 3P) . 



(15) 
(16) 



where = dV/dcj). 



IV. PERTURBATION THEORY IN THE JORDAN FRAME 



In this section, we present the equations that govern the evolution of scalar perturbations in /(i?) theories in the 
Jordan frame. We follow the notation of Kodama and Sasaki [35] for general perturbations and that of Ma and 
Bertschinger [34] in the conformal Newtonian and synchronous gauges. We use Latin and Greek letters to denote 
spatial and 4-vector coordinates, respectively, and the Einstein summation convention is followed. 

As in [35], we may write perturbations of the metric so as to separate the spatial and time dependences. For a 
given wave-number k one can decompose the metric into four time dependent perturbations A, B, Hl and Ht 

goo = -a''il + 2AY), 
goi = -a^BYi, 

9ij = + '^HLY^ij + 2HTYij) , (17) 
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where is the spatial metric and Y = Y(k, x) is the complete set of scalar harmonic functions. Here we consider 
the flat case, for which Y oc exp(i/e • x). The perturbations in the energy momentum tensor are decomposed into 4 
components: density, Sp = p5; velocity, v; isotropic pressure 5P; and anisotropic stress |(p + P)cr, via 



T% = -p[l + SY], 

n = {p + P)iv-B)Yj, 



(18) 



where we use the notation of Ma and Bertschinger [34] for the anisotropic stress (Kodama and Sasaki [35] use the 
anisotropic stress perturbation 11, where PII = P)a). 

In the Jordan frame, matter is minimally coupled and follows the geodesic of the usual metric 17^,,^. Thus, there is 
no dependence of the matter Lagrangian Cm[Xi,gij,v\ on our f{R) modifications, and the conservation equations for 
matter do not differ from the conservation equations of standard general relativity. 



(5' + (1 + w)(kv + 3Hi) + m (^-^ -w]5 = 

w 



SP/5p , 



v' - B') + nil - "iuuMv -B) + {v -B)- —^k^S - k^A + -k^cr = 

I + w I + w 3 



(19) 

(20) 



Two additional independent equations come from the four perturbed field equations. Defining 5R/Y by the following 

expression 



5R 



-6— - SUA' + k^A + kB' + SkHB + mH'j^ + SH'l + 2k'^ IHl + ^] 



we provide all four equations here for completeness, and provide some intermediate results in appendix A. 
— component 



(21) 



QH^A - 2knB - QHH'i^ - 2k^ Hl + 



5R 



Ht 



, 6R 



SR\' 

Y 



- {k^RR + mf^i,) — - SnfRR ( — ) +f'R [GHA -kB- SH'^) = -K^^pS 



(22) 



— i component 



kHA -k{HL + 



1 



5R\' 1 



- Trfc fRR^ +7;kn( Jrr 



SR 

'y 



-M'rA : 



2 2 



ip + P) {v-B) 



(23) 



i — i component 



2(l + /ii) 



A (W^ + 2W') - h^A - \k {B' + KB) - \hB -H'l- 2nH'^ + HA' - (hl + 



+ (27i^ + Ti') Jrr— - /rr (^~y^ " ^/fi-R (^'y' 
+fji (a' + 2HA - IkB - 2H'jA = K^a'^SP 



(24) 



* ~ i 7^ i) component 



-k'^A -k{B'+ KB) +H^+ UH'^ - P i Hl 



U{H^- kB) 



XT) O 

-fcVflK^ - f'R {kB - H'^) = «2a^|(p + P)a 



(25) 
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A. Conformal Newtonian Gauge 



In the conformal Newtonian gauge Ht = B = 0, A = ip and Hl = —(f>, where we have used the notatfon of Ma and 

Bcrtschingcr (Kodama and Sasaki use Hj^ = $ and Ht = '^)- In this gauge, the evolution equations (19) and (20) for 
a cold dark matter (CDM) and radiation ovcrdcnsity may be combined in single second order differential equations, 



6'^ + ns'^ + k^TP-34>" -met)' = 

S'-; + + ^k^^ - 4$" = 



(26) 
(27) 



These, in combination with the following two independent equations, are sufficient to fully specify the evolution 
* ~ i (^7^ i) component 



(l + /fl)(V-0) + /. 



RR 



SR 



-^K^^^i{pi + Pi)(Ti 



— component 



(1 + fn) [2ect> + an i<p' + H^)] + ifRR-H' ^ 



5R 



- {k^fRR + mfnn) — - SHfRR 



SR\' 



where the perturbed expression for the Ricci scalar is now 



+ /'r (6WV + 3(^') = -K^a^i:iPi5i , 



SR 



-6— V - + - 9n<i)' - 30" - 2fc2 
a 



(28) 



(29) 



(30) 



B. Synchronous gauge 

In the sychronous gauge A = B = and, following Ma and Bertschinger, Hl = h/6 and Ht = —3(7? + h/6). To 
completely define the synchronous coordinates we may remove the remaining freedom by specifying that cold dark 
matter particles have zero peculiar velocity v in this gauge. The evolution equations for cold dark matter and radiation 
then reduce to. 



k"^ 4 
'5;' + y<5.-|<5;' = 0. 



(31) 

(32) 



In this gauge it is again possible to use the following two independent equations to fully specify the evolution 
i — j {i ^ j) component 



(1 + /r) [2k^rj - {h" + 6r?") - 2H {h' + 6r?')] - fk {h' + 6r,') - k^fRR^-^ = K^a^SJ:,{pi + Pi)a 
— component 



(33) 



(1 + Jr) [-m' + 2k^ + sjrrH' ^ - {k^fRR + mfRR) ^ - mfRR (^^ j - -fRh' = -^^a^Y^^pA m 



with the perturbed Ricci Scalar now given by 



SR _ 2_ 
Y ~ 



-Hh' + -h" - 2fc2?? 
2 2 ' 



(35) 
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V. STRUCTURE FORMATION IN LATE-TIME DOMINATING f{R) THEORIES 

We now consider in more detail, the solutions to perturbation evolution in scenarios in which the f{R) coupling 
becomes important at late times. Wc first consider evolution in the Einstein frame, since it is intuitively somewhat 
simpler, for pressureless matter, radiation and a scalar field, neglecting the interaction between baryons and radiation 
and then transform results to the Jordan frame using the approach outlined in appendix B. 

As in section III, the Einstein frame variables (except the scalar field) are denoted by a tilde. In the Einstein frame 
we have background equations 

= y + a'V{ct>) + -a'pe + ~a'p,^ (36) 
(j)" + + aV^ = ^Kpa^Pc (37) 

^ p*exp(-^) (38) 

p:^ ^ (39) 

where p*^ is a constant and the scalar field is decomposed into a time dependent background and spatial varying 
perturbation (j){t) + 5(l){x,t). As discussed in [22], a scalar coupling to dark matter does not change the first order 
perturbation equations for CDM if we define the perturbation with respect to p* 

So ^'4 (40) 
00 = ik'v* , (41) 

where 6c is defined consistently with 6c- With the residual freedom in the synchronous gauge we can then set = 
which relates 6c = —^h as in GR. The perturbation in the Jordan frame is then given by 

Sc = 6e+'^nP6c^ (42) 

The perturbation equations in the synchronous gauge are subsequently 

6'^ +n6'^-^n'^{2n^6^ + nc{6c-^Kp5(p) + 2K'^(P'6(t)' - K^V^6(f) = (43) 

6^" + 2H6^' + k^6cl, + a^V,^^6^-cP'6'^- ^n^nc{6c- -KpSij)) = (44) 

~S" + lk% - = 0. (45) 

Radiat/ion dominuLed era: In the radiation era, a a r, the coupling is negligible and the perturbation evolution 
equations become 

S'c + -~Sc -^K ^ (46) 



which yields solutions, fixing 6ci = Sc{Ti) 



K+h%-p: ^ (48) 



6c ^ Sc^[^] (49) 
4 1- 

5-Y ~ ,, ^e- (50) 

^ 3[i + (^] ^ ^ ' 
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For the scalar field the dominant driver is the CDM density fluctuation, and yields, 



2 1 ~ 

-5<t>' + k''5<l> « — f]c4 (51) 

T ZK T 



^ ^ SciT^ {Pj » Pc) (52) 



2,K t'^ Teq 

f3 pi 1 J 

n 



4k pl^ [1 



(53) 



In the radiation era, the coupling is negligible and the Einstein and Jordan frames are comparable. In terms of the 
Jordan frame expansion rate, a, we find 

4 ^ L (-) (54) 



a. 



6-y 4 1 



~5c 3 [1 + (fcrli] 
KS(j) 1 pI 1 fa 



(55) 
(56) 



so that the scalar field perturbation evolves more rapidly than the matter and radiation density contrasts but has a 
smaller initial value. The amplitude of the scalar field perturbation is k dependent, with larger wavenumbers being 
comparatively suppressed. 

At the matter radiation transition the scalar field density fluctuation is approximately 

Because of the suppression of the scalar field fluctuations, the Jordan CDM density fluctuations in the radiation era 
is comparable to the Einstein CDM density fluctuation 

5c(a) = lc + ^Kl36(t) (58) 

=.4.(^)' (59) 

Matter dominated era: In /(i?) models the background evolution approaches an attractor solution in the matter 
dominated era in which, aoct's {a oc rs), which in the Jordan frame corresponds to a oc t^ {aoc t)[24]. In appendix 
C wc highlight how this result is obtained. The scaling regime has a purely kinetic scalar field with = 1/9, and 
f)c = 8/9, leading to 

The matter and scalar field equations are, 



(60) 



K + l-S'c '^^CSc i«:W) + ^-<5</.' = (61) 



f>" + 2m<t>' + eScj) - J- ij^ - M J_ (5^ _ \Kl35cj>) = (62) 

PKT K T'^ 2 



Denoting x = KS4>eq/ PSc,eg, for modes kreq > 1 the scalar field is always going to be negligible to the matter 
perturbation. Likewise for scales for which I/tq < k < l/rgg suppression of the scalar field will set in during the 
matter dominated era. For these instances, where the scalar field is negligible, we find a scaling solution 

Sc{r) = 6,,eq [ — ] (63) 

^C,eq ( ) (64) 
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For modes for which kTeq < 1 i.e. x^q ^1/4 the scalar field approaches equivalence with the dark matter perturbations 
in the matter dominated era. We find a scaling solution of the form 

^-/34 - 25 1 + A(fc,)2 ^^5^ 

For scales with kreg > 1 the scalar field perturbation is suppressed in the matter era and Sc — Sc- For modes kreq < 1, 
the scalar field scaling solution (65) leads to the Jordan frame perturbation being slightly boosted above that of the 
Einstein frame, 

^c= fl + ^ , . A, .o )~^c (66) 



25 1 + ^(fcT) 

this occurs until r ~ 1/fc at which point the scalar field becomes suppressed and 6c — So- 
Accelerating era: During the accelerating era, the Einstein frame CDM density decreases to effectively zero, so that 
if a (X r^/(i+3»"e//)^ where Weff is the effective Einstein frame equation of state 

Sc+7T-4-' T-^c ~ (67) 

5(j)" + 2m(j)' + a'^V,^^5(l) « (68) 

The potential obtained in the Einstein frame typically tends towards an exponential potential at large F ~ 
exp—\K(j). As discussed in the appendix, the accelerating regime is then an attractor in which We// = — 1 + 
and the scalar field perturbation evolves as 5(j) oc with 

6-A2 
P = --. 



2-/JA 




1±1 1-8/ (69) 



For A < v^, p < with a complex amplitude and the scalar perturbations decay yielding 5c — 5c = constant. As a 
result, in the Jordan frame, 5pc — Pc^c decays ~ a~^. In the current era the classical perturbations are still observed; 
however the decay implies that, as the accelerating era continues, the classical perturbations will be smoothed, and 
quantum fluctuations generated in the era will ultimately become important, as the cosmic no hair theorem requires, 
just as in early universe inflation [45]. 

It is common in f{R) theories that the transition from scaling behavior to acceleration occurs significantly earlier 
than in the ACDM scenario. This leads to the late time suppression of the density fluctuations being more pronounced 
with important implications for comparison with observational data. We discuss the evolution and implications of 
structure formation observations for two specific examples below. 



A. Examples: f{R) = -n'^/R and f{R) = -fiiH'i cxp{-R/fi2Hi) 

It is important to note that the results above are true for any form of f{R) where the coupling dictates the expansion 
rate during the matter dominated era. Such behavior is not sensitive to the details of the potential and is expected 
in a large class of f{R), including all explicit functional forms proposed so far in the literature. 

Here, for concreteness, however we provide numerical results for two specific examples, f{R) = —jJL^/R and f{R) = 
~p,iHgcxp{-R/fi2H^), with n"^^' = K^pra{a = l)/m^ = 0.3 and Hq = 70 /kms-'^Mpc-^. 

In figure 1 , we show the evolution for the specific examples in comparison to an equivalent ACDM scenario for two 
different comoving scales k = 10~^ Mpc~^ and k = 10~^ Mpc~^ relevant to galaxy structure and CMB observations 
respectively. One can see the analytical scaling solutions derived above hold in the radiation and matter dominated 
eras. 

The early onset of acceleration in comparison to the ACDM scenario leads to increased relative suppression of the 
large scale density fluctuations and inconsistencies with the galax;y matter power spectrum and ISW eflect [40, 48]. 
In figure 2 we show the effect on the matter power spectrum at 2 = 0. For the same primordial normalization the 
late time evolution leads to an overall large scale suppression of the matter power spectrum in comparison to the 
ACDM scenario. If one arbitrarily renormalizes the spectrum to be in agreement with galaxy matter spectra, then 
the spectrum still shows too great a suppression at large scales to be consistent with the CMB. 
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FIG. 1: Perturbation evolution for [left panel] the f{R) = —fi'^'/R and [right panel] the f{R) = — /ii-ffg exp(— _R//i2-ffo) model 
with = pin/'iHQ — 0.3 and fi'^, and {^11,^2} chosen in each case to give -ffo ~ 70kms~^ Mpc~^ . [Top panels] The density 

fluctuations for the f(R) theory are compared to that for the equivalent ACDM scenario for two comoving scales k — 10~^ Mpc~^ 
(left) and k — 10~"^A/pc~^ (right). The diamond shows the analytic value of Xc,q in the limit of no suppression from kTf,q > 1, 
as described in (57). [Lower panels] The power law evolution of the density fluctuations dS/dlna for the different density 
components. The main scaling solutions are shown by dotted lines, respectively dlnS/dlna = 3 and ciln5/dlna = 2 for the 
scalar field and matter components in the radiation era and dln(5/dlna = 1.5 and dln(5/dlna = for the matter perturbations 
in the matter dominated and accelerated eras. In the figure (E) and (J) denote the Einstein and Jordan frame quantities 
respectively. 



B. Comment on the subhorizon CDM over density in the Jordan frame 

In [40] a potential problem was observed in the conformal Newtonian gauge in the Jordan frame, in that there 
was an apparent zeroing of the CDM matter density fluctuation at critical scales. This is obviously at odds with the 
results of the previous section so we revisit it here to highlight the origins of, and explain, the discrepancy. 

As seen in section IV evolving the perturbation equations in the Jordan frame is significantly more laborious than 
in the Einstein frame. In order to simplify the equations to study the evolution, therefore, in [40] made a quasi-static 
approximation at sub-horizion scales where k/Ti. 3> 1 . In this approximation, for each perturbation variable, one 
assumes that the time derivative is negligible in comparison to its spatial derivative. 

In the subhorizon limit, the matter evolution equation (26) would simplify to 

5" + nS'^ + e^ = . (70) 
If, following [40], we define the parameter Q as 



Q^-2^-^, (71) 
a-' 1 + fR 



the independent components of the Einstein equations (28) and (29) become 
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FIG. 2: The matter power spectrum for ACDM (full black) and [left panel] the f{R) = —fj.'^/R model and [right panel] 
the f{R) — —^iHq exp(— _R//i2-H'o) model for the same normalization (red dashed) and for normalization to give small scale 
agreement between the two models (red dot-dashed) are shown against the SDSS matter power spectrum data [47]. One can 
see that the f(R) model cannot simultaneously give small scale agreement with galaxy matter power spectrum and large scale 
agreement with the CMB. 



where note that n^a^lp^ 
general relativity 



- pc) 7^ 37i^ here. Compare this to the corresponding scale- independent behavior in standard 



(73) 



Inspection of the overdensity evolution equation (72) leads to an identification of a critical value Q = 2/3, at which 
the overdensity Sc is driven to zero with two additional apparent zero density points at Q = 1/2 and Q = 1, in the 
spatial off-diagonal Einstein equations. This apparent zeroing of the CDM overdensity also arises, under the same 
assumptions, in the synchronous gauge. 

However no such zeroing of Sc is seen to arise from the evolution and subsequent conformal transformation of the 
full, unapproximated equations in the Einstein frame. The discrepancy arises because the quasi-static approximation 
used is too aggressive and removes important information about the evolution. We can see this from doing a conformal 
transformation of the Einstein frame CDM perturbation equation (43), for which we obtain. 



6c+n6c-^K^a^{2p-fS^ 



fRR cn\ 2 2/ I N JrR xu 



Ir 



{l + fR)dT V(l + /fl,) 



Irr 



-6R 



Jrr 



(1 + Ir) 



RJr - f 

(l + fR)' 



SR = 



(74) 



The quasi-approximation therefore, in which one neglects time derivatives for /rr, used to infer the zeroing in the 
Jordan frame, is invalid, since zero, first and second time derivatives of fRR/{l + Jr) all come into play in the 
conformally transformed full equation. 
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VI. CONCLUSIONS 



The possibility that cosmic acceleration is our first signal of a far-infrared modification of General Relativity is a 

logical alternative to the cosmological constant and to dark energy models. 

In this paper, we have focused on f{R) modifications to GR, exploring the details of cosmological perturbation 
theory and its implications for the linear theory of structure formation on the universe. The details of such an analysis 
are complicated by the modified dynamics and we find that it is simpler to conduct the analysis in the Einstein frame. 
In particular, we have seen that the use of the 'quasi-static' approximation in the Jordan frame can lead to misleading, 
incorrect conclusions, suggesting that the overdensity is driven to zero at specific scales. This is inconsistent with the 
conclusions obtained when one uses the conformally transformed results from the full equations in the Einstein frame. 

We have shown, considering the evolution in the Einstein and Jordan frames, and general and specific choices 
of gauges, that the evolution of density fluctuations in f{R) gravity leads to predictions that are inc;onsistent with 
cosmological observations. The large scale density fluctuations are suppressed in comparison to small scales, leading 
to an inability to fit both small scale galaxy data and large scale, CMB, data simultaneously. These findings hold 
for the wide class of f{R) that behave like GR at early times and only diverge from GR at low curvatures, in the 
matter dominated era. As such, these results rule out a wide class of modified gravity models aimed at providing an 
alternative explanation to the dark energy problem. 
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APPENDIX A: JORDAN FRAME PERTURBATION EQUATIONS 



We summarize here some of the components that are used in the derivation of the results in section IV. The 
perturbations to the geometric quantities are unmodified in f{R) theories. For the Christoffel symbols we have: 



5r' 



00 



A'Y, Sr^oj =-[kA + HB] Yj , 



Sr'^ij = [-2nA + {k/3)B + 2nHL + H'^] HjY + [-kB + 2nHT + H^] Yij 
SpQO = -[kA + B' + HB] Y\ 

5V'jk = -kHL {S'jYk + 5\Yj - 7,fcr) + UB^jkY' + Ht (Y'j\k + Y\\j - Yju^' 
And for the Ricci scalar and Ricci tensor. 



5R=^ 



-Q—A - mA' + k'^A + kB' + ikUB + mn'^ + m'l + 2k^ [ Hl + ^ 



Y 



SR 



SRki 



- [k'^A - mA' + k{B' + HB) + 2,H'l + mH'j] Y 
-2 +n^^A- HA- + ^A+^{B' + KB) + ^nkB+ 



+H'l + hUH'i^ + 2 



Hl 



SkiY 



+ 
SRoj 



k^A -k{B' + HB) + H'^+ UH't + 2(^+ U^^ Ht - k^ [hl + + n{H'j. - kB) 



(Al) 
(A2) 
(A3) 
(A4) 

(A5) 



(A6) 
(A7) 



— + S - 2knA + 2kH'T + -kH'^ 



Ykj (A8) 
(A9) 



12 



The main additional components that need to be calculated come from the perturbations to the covariant derivative 
terms in (2) 



We are considering a function F = fji where //{ = fit{f^ + Sfji{x, t), 



^(VOVo/fl) 



—^fRR^- - yjRR^- J + ^ (^""^ -kB- 2,Hl 



a 



+ 2^ A - 2AHf'ji + 4^' 



—fRR^-^ - n [fRR^-§^ ' +f'n {^HA -^B- 2nHr. - H', 



k^fRR^S- + I'r {kB - 2HHt - H't) 



Y 



~kf'RH 



Y 



RR 



km 



RR 



SR 



SiiY + 



f'^ikA + HB) 



(AlO) 
(All) 



(A12) 
(A13) 



(A14) 
(A15) 



APPENDIX B: EQUATIONS FOR THE MAPPING OF PERTURBATION THEORY TO THE EINSTEIN 

FRAME 



As reviewed in section III , there is a particular conformal transformation which maps the Jordan frame action into 

a Hilbert-Einstcin action for gravity, with the introduction of a scalar field (f>{x,t). In this Einstein frame, matter 
doesn't fall along the geodesies of the metric, hence the energy-momentum conservation now reads 

V4f^'^+<7)-0 (Bl) 

where T^'^ is the energy-momentum tensor for matter in Einstein frame and Z:^" is the one associated with the scalar 
field. Following the notation of section 111, we use a tilde to indicate quantities in the Einstein frame, (with the 
exception of the scalar field which doesn't have a counterpart in the Jordan frame). 

In a similar fashion, Einstein equations include now terms dependent on the scalar field, as shown in eq.(14). In 
particular, the terms on the RHS of eq.(14) depending on the scalar field, correspond to the energy-momentum tensor 
of (p, and hence we can rewrite eq.(14) as follows 

G^, = 87rG(7;, + i^J (B2) 

Similarly to what done in section IV, we can introduce small perturbations to the metric and the energy-momentum 
tensors, and derive the first order equations by expanding eq.(Bl) and eq.(B2) in the perturbations [38]. 

Alternatively, we can obtain the same equations by mean of the mapping between Jordan and Einstein frame 
described in section III. In this section we will follow the latter method, and derive explicitily the mapping for the 
perturbed quantities from Jordan frame to the Einstein one. Once we have this mapping, it is straightforward to 
derive the perturbed equations by direct substitution into the equations derived in section IV. 

In the Jordan frame, the perturbed metric is written as 

5/.. = 9^°J + K. , (B3) 

where the perturbation h/^v is assumed small. Similarly, in the Einstein frame, the perturbed metric and scalar field 
are written as 

.9^1^ = g'ij^J = 00 + S4)Y , (B4) 

where Scj) is a time dependent small perturbation and Y = Y{x) is the set of harmonic functions. Using (8), we can 
see that at first order 

5^. = gll'J + = e^"^^" + h^. + f3ngl'j6<pY) . (B5) 
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Similarly, given the energy- momentum tensor in Jordan frame, we can use (13) and (8), we can derive the following 
expression for the perturbation to the energy-momentum tensor of matter in Einstein frame 

fi. = fi^io) + st'^ = e-'^f^'"*" (t^^°K + ST^"^ - 2/3kT^(°)^(5(/)F) . (B6) 

Equations (B5) and (B6), together with the background version of (8) and (13), are all is needed to map the perturbed 
expressions in Jordan frame into the ones in Einstein frame, as we show in the rest of this section. 

1. Metric perturbations 

Let's write the perturbed line element in Einstein frame as follows 

ds^ = -a^ (r) (1 + 2AY)dT'^ - {T)BYidTdx' + d^ (r) (7^ + ^HlYih + 2HTYij)dx' dx^ (B7) 

We do not need to consider vector and tensor perturbations as they cannot be generated through a conformal trans- 
formation. 

Using (B5) and (8) we can find the following explicit expressions for the metric elements 

500 = -a^[l + 2AY + PK6(j)] 
90j = -a^BYj 

9ij = a'^bij + C^Hl + pKS(l))^ij + 2HTYij] . 

Comparing the above expressions with (B7) we obtain the prescription to map the metric perturbations from Jordan 
to Einstein frame; specifically: 

A = A- ^kS4> 
B = B 
Hl = Hl- ^i^H 

Hj" = Hx • {■^^) 

2. Matter perturbations 

In the Einstein frame we can write the first order perturbed energy-momentum tensor for matter, analogously to 
how we defined it in the Jordan frame (17), i.e. 

fo = -p[l + 5Y], 

f° = {p + p){v-B)Yj, 

f; = p[S} + ttlYS'j + TTrYj]. (B9) 
Using (B6) and (13) we can find the following explicit expressions for its components 

fo = -p[l + 6Y -2i3Kd(l)Y] 
f? = {p + p){v-B)Yj 
f j = p[S'j + (TTi - 2f3KS(l>)YS'j + TTrYj] . 
Finally, comparing with (B9) we get the following prescription to map the matter perturbations: 

S = 6 + 2(3K6(j) 

V = V 

TTl = TtL + 2^K5(j) 

ITT = TTT ■ (BIO) 

The set of mappings (B8) and (BIO) is all we need to obtain the first order equations in Einstein frame directly from 
the equations in Jordan frame, which we derived explicitely in section IV. 
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APPENDIX C: EINSTEIN FRAME BACKGROUND EVOLUTION SCALING BEHAVIOR 



In this section, for completeness, we demonstrate the Einstein frame scaling behavior that occurs with the onset of 
non- negligible ,f{R) coupling relevant in section V, as pointed out in [42]. Consider the coupled equations (where the 
dot represents d/di) 



— Cf3K^Pr, 


CPnpm 



K 

y 



H = 



+ V + pm+p-y 



+ Pm + 



(CI) 
(C2) 
(C3) 

(C4) 
(C5) 



taking V ~ ^exp(— k/3/U(^) appropriate for many potentials V{(j)) arising out of the conformal transformations of f{R) 
theories for large values of (f). Following [42] we introduce the parameterization. 



K (J) 



K VV K \/ Pj 

^ — ~ ^ , V = — 1=^ , z = — ^ , q; = In a 

The Priedmann equation (C4) allows the matter density to be defined in terms of x, y, z 

1^^ Pm 1 2 2 2 

— — = 1 — X — V — z 

and 



H 



l + x'-y^ + — 



Representing d/da by ', equation (C2) gives 

z' = -z[l-3x^ + 3y'^ - z^] 



In the presence of radiation 



H 

y' = -pyx + y 



= x[j-l^+py^ + C{l-x^-y^-z^) 



If the radiation contribution is negligible one can't use (CIO) and instead one gets 



y = y 



-px + - (1 + - 2/2) 



x' = C-^-Cx^+^-^+y^L-C-^ 



(C6) 

(C7) 

(C8) 
(C9) 

(CIO) 
(Cll) 
(C12) 

(C13) 
(C14) 



Scaling attractors satisfy the constraint x' = y' = z' = Q thereby giving f2^=const. We are interested in the regime 
in which the f{R) coupling (for which C = 1/2) comes to be important in the matter dominated era. From (C13) 
and (C14) we see that the attractor requires y = and x = 1/3 so that il^ = 1/9 and ii;^ = 1 i.e. the scalar field is 
purely kinetic, and We// = 1/9- In this case 



o oc ^5 oc r2 . 



(C15) 
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Considering the Friedmann equation, therefore, one finds a scalar field solution (f) = <j)(j\n.{t/tQ) with t^o = 3//3K. 
Converting to the Jordan frame one obtains 



In the accelerating regime the Einstein frame evolution tends towards and attractor with x = \/^/^ and y = 
\/(l - A2/6) and « so that Weff = 1- X^/S. 
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